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Abstract. We report on the analysis of 100 ks INTEGRAL observations of the Galactic microquasar GRS 1915+105. We focus 
on INTEGRAL Revolution number 48 when the source was found to exhibit a new type of variability as preliminarily reported 
in Hannikainen et al. (2003). The variability pattern, which we name f, is characterized by a pulsing behaviour, consisting of 
a main pulse and a shorter, softer, and smaller amplitude precursor pulse, on a timescale of 5 minutes in the JEM-X 3-35 keV 
lightcurve. We also present simultaneous RXTE data. From a study of the individual RXTEfPCA pulse profiles we find that the 
rising phase is shorter and harder than the declining phase, which is opposite to what has been observed in other otherwise 
similar variability classes in this source. The position in the colour-colour diagram throughout the revolution corresponds to 
State A (Belloni et al. 2000) but not to any previously known variability class. We separated the INTEGRAL data into two 
subsets covering the maxima and minima of the pulses and fitted the resulting two broadband spectra with a hybrid thermal- 
non-thermal Comptonization model. The fits show the source to be in a soft state characterized by a strong disc component 
below ~6 keV and Comptonization by both thermal and non-thermal electrons at higher energies. 
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1. Introduction 

GRS 1915+105 has been extensively observed at all wave- 
lengths ever since its discovery. It was originally detected 
as a hard X-ray source with the WATCH all-sky monitor on 
the GRANAT satellite (Castro-Tirado, Brandt & Lund 1992). 
Apparent superluminal ejections have been observed from 
GRS 1915+105 several times with the Very Large Array in 
1994 and 1995 (Rodriguez & Mirabel 1999) and in 1997 
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with the Multi-Element Radio-Linked Interferometer Network 
(Fender et al. 1999). Both times the true ejection velocity was 
calculated to be > 0.9c. In addition to these superluminal ejec- 
tions at the arcsec scale, GRS 1915+105 has also exhibited a 
compact radio jet resolved at milli-arcsec scales corresponding 
to a length of a few tens of AU (Dhawan, Mirabel & Rodriguez 
2000; Fuchs et al. 2003). Using the Very Large Telescope, 
Greiner et al. (2001) identified the mass-donating star to be of 
spectral type K-M III. The mass of the black hole was deduced 
by Harlaftis & Greiner (2004) to be 14.0 ± 4.4 M Q in a binary 
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orbit with the giant star of 33.5 days. The Rossi X-ray Timing 
Explorer (RXTE) has observed GRS 1915+105 since its launch 
in late 1995 and has shown the source to be highly variable on 
all timescales from milliseconds to months (see e.g. Belloni et 
al. 2000; Morgan, Remillard & Greiner 1997). Belloni et al. 
(2000) categorized the variability into twelve distinct classes 
which they labelled with Greek letters (two of those repre- 
sented by semi -regular pulsing behaviour), and identified three 
distinct X-ray states: two softer states, A and B, and a harder 
state, C. The source has been detected up to ~ 700 keV dur- 
ing OSSE observations (Zdziarski et al. 2001). For a recent 
comprehensive review of GRS 1 9 1 5 + 1 05 see Fender & Belloni 
(2004). 

The European Space Agency's INTErnational Gamma-Ray 
Astrophysical Laboratory (INTEGRAL) is aimed at observing 
the sky between ~ 3 keV and 10 MeV (Winkler et al. 2003). 
The INTEGRAL payload consists of two gamma-ray instru- 
ments - the Imager on Board the INTEGRAL Spacecraft, IBIS 
(Ubertini et al. 2003), and the Spectrometer on INTEGRAL, 
SPI (Vedrenne et al. 2003); two identical X-ray monitors - the 
Joint European X-ray monitor, JEM-X (Lund et al. 2003); and 
an optical monitor (Mas-Hesse et al. 2003). 

GRS 1915+105 is being observed with INTEGRAL as part 
of the Core Program and also within the framework of an 
Open Time monitoring campaign. Here we present the results 
of the first of our AO-1 Open Time observations which took 
place during INTEGRAL Revolution number 48. The source 
was found to exhibit a new type of variability not documented 
in Belloni et al. (2000). Preliminary results were described 
elsewhere (Hannikainen et al. 2003). We also have a simul- 
taneous observing campaign with RXTE to concentrate on tim- 
ing analysis (Rodriguez et al. 2004a). The main goal of this 
INTEGRAL/RXTE program was to study the disc-jet connec- 
tion in GRS 1915+105 and to study the variability classes in 
the hard X-ray regime. As the source behaviour is very un- 
predictable, we were not able to study the disc-jet connec- 
tion in this particular case, nor any of the already known vari- 
ability classes. Instead, (as we show below), the source ex- 
hibited a new type of variability and we discovered a new 
source (IGR J19140+0951) during the observation. In this pa- 
per, we conduct a more detailed analysis of the high energy be- 
haviour of GRS 1915+105 during Revolution 48. In Sectional 
we discuss the observation and data reduction methods, while 
in Section [3] we re-introduce the new type of variability. In 
Section|4]we present the RXTE results, while in Section[5]we 
discuss the spectral analysis of the INTEGRAL data. Finally, in 
Section[6]we present a summary of the paper. 

2. Observations and data reduction 

As part of a monitoring programme which consists of six 100- 
ks observations, INTEGRAL observed GRS 1915+105 for the 
first time for this campaign on 2003 Mar 6 beginning at UT 
03:22:33 during INTEGRAL'S Revolution 48 for 100 ks. The 
INTEGRAL observations were undertaken using the hexago- 
nal dither pattern (Courvoisier et al. 2003), which consists of 
seven pointings around a nominal target location (1 source 
on-axis pointing, 6 off-source pointings, each 2° apart, with 
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Fig. 1. The RXTE/ ASM one-day average 1.3-12 keV lightcurve 
showing the epoch of the INTEGRAL observation and the si- 
multaneous pointed RXTE observation (dashed line). 



each science window of 2200 s duration). This means that 
GRS 1915+105 was always in the field-of-view of all three 
X-ray and gamma-ray instruments (JEM-X, IBIS and SPI) 
throughout the whole observation. 

Figure G] shows the RXTE/ASM one-day average 1.3- 
12 keV lightcurve. The vertical line indicates the date of the 
INTEGRAL and the simultaneous RXTE pointed observations. 

2.1. INTEGRAL data reduction 

All data were reduced using the Off-line Scientific Analysis 
(OSA) version 4.0 software. The JEM-X data were reduced fol- 
lowing the standard procedure described in the JEM-X cook- 
book. Based on the JEM-X 3-35 keV lightcurve first presented 
in Hannikainen et al. (2003) and reproduced here (Fig.|2j, we 
defined Good Time Intervals (GTIs) for the spectral extrac- 
tion process. We created two GTI files: one for data below 
70 counts s _1 (which corresponds to ~ 0.5 Crab in that energy 
range) and one for data above 70 counts s _1 . We then ran the 
analysis using the two GTI files from which we extracted two 
total spectra. For the modelling of the spectra, we ignored data 
above 30 keV for the spectrum of the bright parts and above 
25 keV for the faint parts. We assumed systematic uncertain- 
ties as shown in Table [2 (P. Kretschmar & S. Martinez Nunez, 
priv. comm.). 

Concerning IBIS, we reduced the data from the top layer 
plane detector, ISGRI (Lebrun et al. 2003). We made a run of 
the software up to the IMA level, i.e. the production of images. 
During the run we extracted images from individual science 
windows in two energy ranges (20^40 keV and 40-80 keV) as 
well as a mosaic in the same energy ranges. As with JEM-X, 
the GTIs were given as input to OSA to extract the ISGRI spec- 
tral products. Rather than using the standard spectral extrac- 
tion, we then re-ran the analysis software up to the IMA level 
in ten energy bands and then estimated the spectra. The source 
count rate and error were then extracted from each mosaic at 
the position of the source, from the intensity and variance map 
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Fig. 2. The JEM-X 3-35 keV lightcurve for the whole of Revolution 48. The time bin is 8 s. As can be seen, the source varied 
between ~ 0.25-2 Crab. Reproduced from Hannikainen et al. (2003). Shown on the plot is our simultaneous RXTE coverage. 
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Table 1. Level of systematic uncertainty applied to the spectral 
channels of the JEM-X spectra, and energy channel correspon- 
dence. 



respectively (see the Appendix in Rodriguez et al. 2004b for an 
in-depth description of this method). We further assumed 8% 
systematics in each channel (e.g. Cadolle-Bel et al. 2004). 



To reduce SPI data a pre-release version of OSA 4.0 was 
used, which incorporates two different analysis methods: The 
default pipeline is optimised for speed, but cannot handle user 
defined GTIs (Knodlseder et al. 2004). Hence the alternative 
pipeline, which is more closely based on the OSA 3 release, 
was adopted but using the more recent calibration data from 
OSA 4.0. The background flux was estimated using the Mean 
Count Modulation (MCM) method and spectra were extracted 
from the reduced data, using the same input catalogue as the 
ISGRI analysis, with SPIROS v8.0.1 (Skinner & McConnell 
2003). 

The combined JEM-X, ISGRI and SPI spectra resulted in 
two broadband spectra: one for the high part of the JEM-X 
lightcurve, which we will refer to as pulse maxima and one 
for the low part, corresponding to the combined pulse min- 
ima. These two spectra are used for the broadband modelling 
in Section|5] 

2.2. RXTE data reduction 

During the last third of the INTEGRAL observations we ob- 
tained simultaneous RXTE observations (~ 32 ks). We re- 
duced and analysed the data from the RXTE Proportional 
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Fig. 3. a The RXTEfPCA 2.5-15 keV lightcurve (with time bin 
of 1 second) extracted from all PCU and layers. The observa- 
tions start on MJD 52704.99. Marked in the plot are the times 
of the pulses dealt with in Figs. |4] |5] |6] |8] and in Section [4] 
Each segment in the lightcurve represents one RXTE orbit, b 
An expanded plot of the pulses. 

Counter Array using the LHEASOFT package v5.3, follow- 
ing the standard steps explained in the PCA cookbook; see 
Rodriguez, Corbel & Tomsick (2003) for the details of the se- 
lection criteria. We first extracted a 1-s resolution lightcurve 
between 2.5-15 keV from the binned data, from all available 
Proportional Counter Units (PCU) turned on during individ- 
ual ~ 90 min RXTE revolutions. The background-corrected 
lightcurve is shown in Fig. [3] Since we were interested in study- 
ing the spectral variations of the source during this observation, 
we took advantage of the high sensitivity of the PCA to ex- 
tract spectra on the smallest timescale allowed by the standard 
2 data. In a manner similar to a previous work (Rodriguez et 
al. 2002), we extracted PCA spectra and background spectra 
every 16 s from the top layer of PCU and 2. The responses 
were generated with pcarsp (V10.1). The resulting 16-s RXTE 
spectra are used to study fast spectral variability in Section|4] 

3. A new type of variability? 

Fig.0shows the entire JEM-X 100 ks lightcurve. Throughout 
most of the lightcurve (the high fast flaring parts, e.g. from 
MJD 52704.15-52704.41) we see a novel class of variabil- 
ity, characterized by 5-minute pulses, not observed before 
(Hannikainen et al. 2003). The entire JEM-X lightcurve is dom- 
inated by these 5-minute pulses, as evidenced by the 3 mHz 
quasi-periodic oscillation which resulted from a Fourier trans- 
form the whole JEM-X lightcurve (Hannikainen et al. 2003). 
Although this kind of pulsed variability resembles the p- 
heartbeat and k oscillations described in Belloni et al. (2000), 
the oscillations presented there are more uniform and occur on 
shorter timescales. In Hannikainen et al. (2003), we also pro- 
duced a colour-colour diagram for the RXTE data in the same 
manner as in Belloni et al. (2000) and showed that the data 
points fell in a different part of the colour-colour diagram as 
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Fig. 4. The JEM-X colour-intensity diagram from the whole 
100 ks observation. The colour was constructed from the 10- 
35 keV and 3-10 keV lightcurves. The maxima of the flares are 
in the right-hand edge of the plot, while the minima are in the 
left-hand edge. This confirms that the rising phases are harder 
than the declining phases. Each point is one time bin (10 s) 
and represents the averaged mean value from the 300 JEM-X 
flares; the cross shows the average standard deviation for the 
mean value of each individual point throughout the loop - this 
includes both the variability within each point and the error (i.e. 
it is not an error bar). 



compared to the other classes. We call this new type of vari- 
ability 

To study the behaviour of the source in this new variability 
class, we produced a colour-intensity diagram from the whole 
100 ks JEM-X lightcurve (Fig.©. The JEM-X colour was con- 
structed from the 3-10 keV and 10-35 keV lightcurves. Each 
point is one time bin of 10 s and represents the averaged mean 
value from 300 JEM-X pulses. Points from the pulse-maxima 
are situated in the right-hand side of the plot, while points from 
the minima are on the left-hand side. Indicated in the figure are 
the rising and declining phases of the pulses. The cross shows 
the standard deviation of the mean value of each point - this 
includes both the real variability and the error, and hence is not 
an error bar in the conventional sense. The figure shows not a 
simple flux-hardness correlation, but rather some kind of limit 
cycle/hysteretic behaviour with the rising phases being harder 
than the declining phases. This is opposite to what has been ob- 
served in other similar classes (e.g. class p). In the next section, 
using the RXTE data, we take a more detailed look at this be- 
haviour, starting with an examination of the individual pulses. 

4. Pulses and variability in the RXTE data 

During some of the RXTE pointings, the 5-minute pulses are 
clearly detected also by RXTE/PCA (Fig.fJk). In one segment 
in particular (Fig. |3j)) we see nine consecutive pulses, repre- 
sentative of the variability pattern in the majority of the JEM-X 
lightcurve. Because of the better temporal resolution and higher 
collecting area of RXTE/PCA compared to INTEGRAL/JEM- 
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Fig. 5. The mean of nine RXTE/PCA pulses - the times covered 
by these nine pulses is shown in Figs. [3] and The top (red) 
line shows the 2-5 keV lightcurve, the middle (cyan) the 5- 
13 keV and the bottom (blue) the 13^40 keV lighcurve. The ris- 
ing phases are shorter and harder than the declining phases. The 
latter eight pulses were shifted onto the first pulse as described 
in the text, using the maximum of the 2-5 keV lightcurve as the 
reference point. 

X, we can study the variability of the source on a timescale of 
1 s with good enough statistics. This time bin is appropriate be- 
cause Belloni et al. (2000) used a time bin of 1 s to make their 
colour-colour diagrams, hence allowing for a direct comparison 
between this work and theirs. Another reason why a 1 -s binning 
was chosen is that it is not of too high resolution, and therefore 
does not eliminate too many points which would render the plot 
unclear, but at the same time is not so low so as to smear out any 
variability present. The RXTE/PCA lightcurve was divided into 
three energy bands (2-5, 5-13, and 13^40 keV) and smoothed 
with a boxcar average of 5 time bins. Figure[5]shows the mean 
of the nine RXTE/PCA pulses from Fig. 3b. The pulses con- 
sist of the main pulse (with the rising phase being shorter and 
harder than the declining phase) and a precursor pulse, which 
is shorter, softer and of smaller amplitude than the main pulse. 

The times of the maxima of the 2-5 keV pulses were de- 
termined from the smoothed lightcurve. The eight latter pulses 
were shifted onto the first pulse using the times of the maxima 
as the reference. This was repeated for the other two energy 
bands, using the times of the maxima of the 2-5 keV pulses. 
From the resulting average pulse shape, we can confirm that 
the rising phase of the (main) pulse is harder than the declining 
phase. In addition, the rise is also shorter than the decline. In 
each of these nine pulses, there is also a precursor seen only in 
the softer X-ray bands. 

To further study the indicated limit cycle, we constructed 
a colour-colour diagram from the RXTE/PCA data, shown in 
Figure [6] The X-ray colours were defined as HR1 = B/A and 
HR2 = C/(A + B), where A, B, and C are the counts in the 2-5 
keV, 5-13 keV and 13-40 keV, respectively, following Belloni 
et al. (1997). In Fig. [6] the dots represent the whole RXTE ob- 
servation (including the non-pulsed part), while the solid line 
shows the cycle traced by the pulses (Fig.^). The main aim 
of the figure is to show the small range of variability of the 5- 
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Fig. 6. Colour-colour diagram showing the the RXTE/PCA 
from this study. The dots represent the whole observation, 
while the solid line superposed is from the pulses (Fig.[5J>). The 
curved dotted lines correspond to the colours expected from 
a Comptonizing spherical corona with a given radial optical 
depth T50 (i.e. optical depth for the assumed electron tempera- 
ture of 50 ke V) and varying values of the inner disc temperature 
T m . The solid lines show the colours expected for the model for 
given T m and varying t 5 q. The pulse cycle goes clockwise. 

minute pulses. The pulse cycle goes clockwise as indicated by 
the arrows. The maxima of the pulses are at the upper-right of 
the pulse cycle curve, while the minima are at the lower left. 

Superposed in Fig. [6] are model curves from a study by 
Vilhu & Nevalainen (1998), calculated as explained below. 
Although a direct comparison of the parameters is not possible 
due to the RXTE/PCA gain evolution through the years, 
resulting in a change of the energy channel correspondence, 
we simply want to illustrate qualitatively the position in the 
colour-colour diagram of this present observation. In order 
to have an idea of what physical parameters the observed 
colours correspond to, we simulated the colours predicted by a 
Comptonization model. We assumed the Comptonizing source 
to be a sphere surrounded by a cooler disc with the inner radius 
equal to that of the sphere (Poutanen, Krolik, & Ryde 1997). 
Seed photons coming from the disc are characterized by the 
inner disc temperature T m . The electrons in the cloud were 
assumed to be thermal with temperature kT = 50 keV. The 
Thomson optical depth of the cloud along the radius is T50, 
with subscript 50 meaning that this r corresponds to the 50 keV 
electrons. For a different T, similar Comptonization spectra 
are produced by a cloud of optical depth t = T50 x 50keV/fcr, 
because the slope of the Comptonization spectrum, T, is 
a function of the Kompaneets parameter, y = ArkT /m e c 2 , 
and r = \y- 219 (Beloborodov 1999). For the given T m and 
T50, we compute the Comptonization spectrum using the 
method of Poutanen & Svensson (1996). Then we obtain 
the colours expected from the model spectrum using the 
responses of RXTE/PCA and plot them on the colour-colour 
diagram (Fig.|6j. We see that the pulses discussed in this paper 
correspond to optical depth r 0.7 and inner disc temperature 
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of T m ~ 1.3-1.4 keV. This temperature differs slightly from 
that obtained using a black body + power law model, because 
the Comptonization spectrum is not a power law close to the 
seed photon energies (see e.g. Vilhu et al. 2001). 



4.1. Hysteresis or limit-cycle oscillations? 

The colour-intensity diagram in Fig. |4] and the colour-colour 
diagram in Fig. |6] both show loop-like features, which could 
be due to a hysteretic process or limit-cycle oscillations. The 
hysteresis phenomenon has been suggested to explain the loop 
features in the X-ray colour-colour diagram of black hole low- 
mass X-ray binaries (Miyamoto et al. 1995). One recent ex- 
ample is GX 339-4 (see Zdziarski et al. 2004). Hysteresis 
can occur on a large range of timescales, and for accreting 
black holes the timescale is attributed to the finite time re- 
quired in the building up and dissipation of a viscous, optically 
thick accretion disc. The accretion disc gives rise to low-energy 
photons, which act as coolants for the hot plasma via inverse 
Comptonization. During the accretion disc build-up process, 
the source luminosity increases but the X-ray spectrum soft- 
ens. On the other hand, when the accretion disc dissipates (e.g. 
due to evaporation), the source luminosity decreases and the 
X-ray spectrum hardens. The process is expected to occur in 
an extended region at a substantial distance from the central 
black hole. The duration of the two hysteretic states and their 
transition times are long, and for low-mass X-ray binaries the 
timescale would be a year or more. 

In GRS 1915+105, the timescale for developing the loop 
features is relatively short - roughly hundreds of seconds. The 
observation is not consistent with the scenario requiring an 
extended accretion disc, the building up and the recession of 
which drive the hysteretic cycle. One may argue that the pro- 
cess could occur in a smaller inner disc region, and this nat- 
urally shortens the corresponding timescales. When we exam- 
ined the loop in Fig. |4] we find spectral hardening accompa- 
nied by an increase in the X-ray photon counts, contrary to the 
properties of the hysteresis model for the GX 339-4 behaviour. 
We therefore conclude that the processes behind the loops in 
the colour-colour diagrams of GRS 1915+105 and GX 339-4 
are not identical. As an alternative, we interpret the loops in 
the GRS 1915+105 data as limit-cycle oscillations, which are 
probably due to delayed feedback in the inner accretion disc. 
A possible explanation for the simultaneous rise in the photon 
counts and spectral hardening is an increase in the size of the 
hard X-ray emission region which produces the Comptonized 
hard X-rays. This region could be a hot disc corona or some 
non-thermal plasmas in the vicinity of the accreting black hole. 

4.2. Fast spectral variability 

To track variability on short timescales, the 16-s PC A spectra 
were fitted individually in XSPEC VI 1.3 (Arnaud 1996) with 
a simple model of absorbed power law plus multicolour disc 
blackbody (MCD, Mitsuda et al. 1984), and a Gaussian. The 
absorption column density Nh was fixed at 5 x 10 22 cirT 2 . The 



power law photon index F, inner disc temperature kT m and in- 
ner disc radius i?; n are plotted in Fig.0 in addition to the actual 
lightcurve in flux units. All bad fits (points with x 1 > 50 for 
~ 41 d.o.f.) were removed and the data were smoothed with 
a boxcar average of 5 time bins to highlight the general trend 
in the data sets. We clearly see amplitude variations of the fit 
parameters related to variations in the 2-20 keV flux. The in- 
ner disc radius is obtained from the normalisation, N, of the 
MCD model following R m - d\o x y/N/ cos 6, where 6 is the 
inclination angle, and d\o is the distance in units of 10 kpc. For 
plotting purposes, we used 9 = 70° and dio = 1 (Chapuis & 
Corbel 2004 quote a distance range of 6-12 kpc). No colour 
correction (e.g. Shimura & Takahara 1995) has been applied 
to those values. We remark that sometimes the inner disc tem- 
perature shows high values {kT m > 1 .5 keV), along with very 
low values of the inner disc radius. Those very high kT m spikes 
are indicative of times when the MCD model is inappropriate 
(Merloni, Fabian & Ross 2000). However, the simple model 
we used in our fitting allows us to track the general spectral 
variations of the source on these short timescales, and even if 
the MCD parameters are unphysical during the spikes, they are 
indicative of underlying changes in the source spectrum. 

Fig- ID shows correlation plots between F, T\ n , R[ n and the 
20-50 keV flux against the 2-20 keV flux for the whole set of 
RXTE/PCA observations (a-d) and only for the subset of pulses 
(e-h) fitted with the simple MCD+power law model. The over- 
all RXTE/PCA data sets (panels a-d) represent a mixture of 
source behaviour (see the lightcurve in Fig. [3] for example), 
as the RXTE/PCA observations covered times when the source 
behaviour did not necessarily exhibit the new variability pat- 
tern. Fig.[8k shows that there is no obvious correlation between 
the photon index and the 2-20 keV flux throughout the whole 
observation. Panel b shows that there is the tendency for T m 
to increase with increasing flux, while panel c shows that /?; n 
decreases with increasing flux. The data sets from the pulses 
are slightly different from the overall data set which includes 
the non-pulsed part. To begin with, Fig. [3^ shows that on the 
whole the photon index decreases with increasing flux, which 
supports the scenario presented, for example, in Fig.[3]that the 
maxima are harder than the minima. Fig. |8f shows that there 
is no substantial variation in T m with flux, indicating that the 
pulses are largely isothermal. Again, Fig. [8J* shows no obvi- 
ous correlation between 7?; n and flux, i.e. Ri n does vary, but not 
coherently with the flux. 

On the other hand, Fig.|SJl and h show a very strong linear 
correlation between the flux in the 2-20 ke V range and the 20- 
50 keV range, both for the overall data set and for the subset of 
pulses. Taken together with the photon index-flux anticorrela- 
tions, this implies that there is no pivoting in the RXTE energy 
range, but the flux at higher energies increases more than at 
lower energies. 

5. Broadband spectral modelling 

5.1. Spectral Model 

In order to conduct a more detailed analysis of the spectrum, we 
used the co-added spectra, extracted from the high and low flux 
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Fig. 7. The results from the RXTE spectral analysis. Each column represents one RXTE orbit - the fourth column from the left is 
from the time of the pulses examined in more detail. The top row reproduces the PCA 2-20 keV flux. The second row from the 
top shows the evolution of the photon index F throughout the observations, the third row shows the inner disc temperature (keV), 
while the bottom row shows the inner disc radius (km). All results are from the fits to the 16-s spectra. Bad fits, i.e. data points 
with^ 2 > 50 (for ~ 41 d.o.f.), were removed. The remaining data and fit results were smoothed with a boxcar average of 5 time 
bins to show the general trend throughout the data sets. 



parts in the INTEGRAL lightcurve as described in Section 13 
This provides us with two broadband spectra (JEM-X + ISGRI 
+ SPI) in the energy range of 3-200 keV, for the pulse max- 
ima and minima separately. The two spectra were fitted us- 
ing XSPEC VI 1.3 (Arnaud 1996). We interpret the spectra 
in terms of Comptonization of soft X-ray seed photons, as- 
sumed here to be a disc blackbody with a maximum temper- 
ature, Tin. We use a Comptonization model by Coppi (1992, 
1999), eqpair. This model has previously been used to fit X- 
ray spectra of GRS 1915+105 (based on /ttTE/PCA/HEXTE 
and CGRO/OSSE data) by Zdziarski et al. (2001), Cyg X-l 
(Poutanen & Coppi 1998; Gierlihski et al. 1999), as well as 
INTEGRAL and RXTE data of Cygnus X-3 (Vilhu et al. 2003). 

The electron distribution in the eqpair model can be 
purely thermal or hybrid, i.e., Maxwellian at low energies 
and non-thermal at high energies, if an acceleration process is 
present. The electron distribution, including the temperature, 
T e , is calculated self-consistently from the assumed form of the 
acceleration (if present) and from the luminosities correspond- 



ing to the plasma heating rate, Lh, and to the seed photons irra- 
diating the cloud, L s . The total plasma optical depth, t, includes 
contributions from electrons formed by ionization of the atoms 
in the plasma, t; (which is a free parameter) and a contribu- 
tion from e* pairs, r — Tj (which is calculated self-consistently 
by the model). The importance of pairs and the relative im- 
portance of Compton and Coulomb scattering depend on the 
ratio of the luminosity, I, to the characteristic size, r, which 
is usually expressed in dimensionless form as the compact- 
ness parameter, { = Lcrj /(rm e c 3 ), where cr T is the Thomson 
cross section and m e is the electron mass. The compactnesses 
corresponding to the electron acceleration at a power law rate 
with an index, Fi n j and to a direct heating (i.e., in addition to 
Coulomb energy exchange with non-thermal e ± and Compton 
heating) of the thermal e ± are denoted as ^ nt h and (±, respec- 
tively, and = ^ nt h + Ah- Details of the model are given in 
Gierlinski et al. (1999). 

Following Zdziarski et al. (2001), we assume here a con- 
stant {$ = 100, compatible with the high luminosity of GRS 
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Fig. 8. Panels a -d show T, T m , R m and the 20-50 keV flux against the 2-20 keV flux from the whole set of the RXTE/PCA 
observations, while panels e-h show the same but only for the pulses. 



Table 2. Model parameters 3 for the best-fit models to the broadband spectra for the co-added pulse maxima and minima. 
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0.50 
0.34 
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1 5 +0 - 15 

i- -0.15 

2.0_ 10 


2.2 x 10- 3 

1.3 x 10- 2 


5.8 x 10~ 8 

2.9 x 10- 8 
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a The uncertainties are for 90% confidence, i.e., Ay 2 = 2.71. 

b Calculated from the energy and pair balance, i.e., not a free fit parameter. 

c Assumed < 2 in the fits. 

d The bolometric flux of the unabsorbed model spectrum. 



1915+105. For example, for 1/2 of the Eddington luminosity, 
Le, and spherical geometry, the size of the plasma corresponds 
then to r ~ 100GM/c 2 . We note, however, that the dependence 
of the fit on 4 is rather weak, as this parameter is important 
only for e ± pair production and Coulomb scattering, with the 
former not constrained by our data and the latter important only 
at ^ < 1 or so (Gierlihski et al. 1999). 

We include Compton reflection (Magdziarz & Zdziarski 
1995), parametrized by an effective solid angle subtended by 
the reflector as seen from the hot plasma, £1, and an Fe Ka flu- 
orescent line from an accretion disc assumed to extend down 
to 6GM/c 2 . We assume a temperature of the reflecting medium 
of 10 6 K, and allow it to be ionized, using the ionization calcu- 
lations of Done et al. (1992). We define the ionizing parameter 



as ^ on = 47rF lon /n, where Fj on is the ionizing flux and n is the 
reflector density. As the data poorly constrain £j on , but clearly 
require the reflector to be ionized, we freeze it to a value of 100 
erg cm s~', in the middle of the confidence intervals for both 
fits. 



We further assume the elemental abundances of Anders & 
Ebihara (1982), an absorbing column of A^h > 1.8 x 10 22 crrT 2 , 
which is an estimate of the Galactic column density in the di- 
rection of the source by Dickey & Lockman (1990), and an 
inclination of 66° (Fender etal. 1999). 
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5.2. Results 

The two spectra are shown in Fig. [5] together with the best- 
fit model. The parameters for the best-fit models are shown in 
Table|3 and Fig.llOlshows the spectral components of the fits 
for the two spectra. Note that the small values for^ 2 are due to 
the systematics assumed for the INTEGRAL data. 

Figure [9] shows that the main difference between the pulse 
maxima and minima is in the observed flux, which is higher 
by a factor of ^2 during maxima. Both spectra show that GRS 
1915+105 during Revolution 48 was in a soft state, with spec- 
tra similar to that of state A in Sobolewska & Zycki (2003) 
and in between the softest spectrum in the variability class C/x 
and that in the B/y class, as observed by RXTE and CGRO in 
Zdziarski et al. (2001). 

The fits imply (Fig. llOt that during pulse maxima as well 
as minima the spectrum at 56 keV is dominated by the unscat- 
tered disc emission, with Comptonized emission dominating 
only at higher energies, and including a significant contribu- 
tion from non-thermal electrons. The figure shows scattering 
by the thermal and the non-thermal electrons separately. This 
is obtained by first calculating the scattering spectrum by the 
total hybrid distribution, consisting of a Maxwellian and a non- 
thermal high-energy tail. Then, the non-thermal compactness 
is set to zero, ( nt ^ = 0, leaving us with only thermal plasma. 
However, since the plasma parameters, r and kT e , are deter- 
mined by the model self-consistently, they change after setting 
^nth = 0. Thus we have to adjust the model input parameters 
so as to recover the original r and kT e . Since there is substan- 
tial pair production in the hybrid model but virtually none in 
the thermal model, we can set the (input) parameter r; equal to 
the (output) r of the hybrid model. We then need to adjust the 
plasma cooling rate, determined by (h/€ s . Thus, we adjust that 
ratio to a value at which kT e is equal to the original value of the 
hybrid model. At this point, we calculate the thermal scattering 
spectrum, and subtract it from the total scattering spectrum to 
obtain the non-thermal component. A Compton reflection com- 
ponent including a strong Fe Ka line is also required by the fits 
to both spectra. 

We see, however, that there is significant softening of the 
Comptonized part of the spectrum during pulse minima. This 
is reflected in the ratio, being significantly smaller. The 

ratio (h/{s, or, equivalently, L^/L s , is between the power sup- 
plied to the electrons in the Comptonizing plasma and that in 
soft disc blackbody photons irradiating the plasma. Our fits do 
not require an additional component in soft photons not irradi- 
ating the hot plasma and thus L s corresponds to the luminosity 
in the entire disc blackbody emission. 

As implied by the value of the bolometric flux in Table [2] 
(see also Fig. [9ji, the value of L s (which dominates the bolo- 
metric flux) decreases by a factor of ^ 2 during minima. Given 
the corresponding decrease of Lh/L s by a factor of ^ 2, we find 
that Lh decreased by ~ 4. Conversely, there is an increase in 
the relative power supplied to the electrons in the plasma dur- 
ing pulse maxima by a factor of ~2. 

It should be noted that the spectrum from the pulse maxima 
contains both the end of the rising phase and the beginning of 
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Fig. 9. Deconvolved spectra of GRS 1915+105 during 
Revolution 48 from the pulse maxima (top) and minima (bot- 
tom) separately. The cyan, magenta and blue datapoints are 
from the JEM-X, ISGRI and SPI, respectively. The ISGRI and 
SPI data were renormalized to the JEM-X. The dashed and 
solid curves represent the best-fit models to the observed and 
unabsorbed spectra respectively, with the parameters given in 
Table|3 
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Fig. 10. Spectral components of the fits to the pulse maxima 
(upper panel) and minima (lower panel). The dotted and long- 
dashed curves show the unscattered blackbody, and Compton 
scattering from thermal (red) and non-thermal electrons (cyan), 
including a component from e ± pair annihilation, important 
around 511 keV. The short-dashed (green) curve shows the 
component from Compton reflection, including the Fe Ka line. 
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the declining phase of each pulse, and thus it is an average of 
the harder rising and the softer declining phases. 

The other best-fit parameters for the two spectra are consis- 
tent with being the same within the measurement uncertainties. 
Given the importance of e ± pair production, the contribution 
to the total optical depth, r, from the ionization electrons is 
relatively weakly constrained, especially in the high-flux state, 
when it is consistent with being null (i.e., with all scattering 
plasma formed by pair production). The corresponding range 
of the total optical depth, t, is much narrower; however, the 
eqpair model does not allow the uncertainties in r to be cal- 
culated directly as it is not a free parameter. 

6. Summary 

We have presented INTEGRAL and RXTE data for 
GRS 1915+105 from INTEGRAL'S Revolution 48 in 2003 
March. The source was found to exhibit a novel type of vari- 
ability which we call The continuous 100 ks observations, 
due to INTEGRAL'S eccentric orbit, allowed us to classify this 
new variability that dominated the entire JEM-X 3-35 keV 
lightcurve, and is characterized by repeated 5-minute pulses. 
The simultaneous RXTE/PCA observations showed in one 
pointing a series of nine ~ 5 minute pulses, corresponding to 
the pattern seen by JEM-X. By examining the source behaviour 
in terms of flux-hardness correlations as well as individual 
pulse shape, we have shown that the rising phases of the pulses 
are harder and shorter than the declining phases, contrary to 
what is observed in the p and k classes of Belloni et al. (2000) 
which this new variability class otherwise resembles. Simple 
modelling of the fast spectral variability confirms a negative 
correlation between hardness and flux but shows no obvious 
correlation between flux and inner disc temperature or radius, 
suggesting that the pulses are largely isothermal. Loop-like 
features, traced out by the pulses in the colour-colour diagram 
are interpreted as limit-cycle oscillations probably due to 
delayed feedback in the inner accretion disc. 

Broadband spectral modelling of co-added spectra for the 
pulse maxima and minima respectively, using eqpair, showed 
that the spectra of maxima and minima both belong to spectral 
state A, according to the classification by Belloni et al. (2000), 
and were essentially similar, except for a significant softening 
of the Comptonized part of the spectrum during minima. The 
similarity of the spectra was also confirmed in the plot of the 
2-20 keV against the 20-50 keV PCA fluxes which was essen- 
tially linear, both throughout the whole observation and solely 
during the pulses. 

We have also seen that there is essentially no correlation 
between T m and the flux during the pulses, again supporting 
the suggestion that they are isothermal. This implies that this 
new variability pattern in itself is largely isothermal by nature. 
The fact that the spectra do not change between the pulse max- 
ima and the pulse minima (except for the flux level) and that 
the pulses are described by a fast rise and a slow decline in- 
dicates that the variability is not representative of oscillations 
between spectral states (as e.g. class p in Belloni et al. (2000)). 
We attributed the variability to limit-cycle oscillations which 
are probably due to delayed feedback in the inner accretion 



disc. However, the variability could also result from flares in 
the accretion disc arising from instabilities in the accretion rate. 

Recently, Tagger et al. (2004) have shown, based on the 
work by Fitzgibbon et al. (priv. comm.), that GRS 1915+105 
always seems to follow the same pattern of transition through 
all 12 classes as defined by Belloni et al. (2000). This may sug- 
gest that GRS 1915+105 simply displays a repeating contin- 
uum of all possible configurations and that we caught it in an 
intermediate class between two previously known ones. This 
implies that more variability classes remain to be observed and 
that, with continued monitoring of GRS 1915+105, we should 
be able to observe other classes of variability behaviour and 
gain a better understanding of this mystifying source. 
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